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Recent type la supernovas data seemingly favor a dark energy model whose equa- 
tion of state w{z) crosses —1 very recently, which is a much more amazing problem 
than the acceleration of the universe. In this paper we show that it is possible to 
realize such a crossing without introducing any phantom component in a Gauss- 
Bonnet brane world with induced gravity, where a four dimensional curvature scalar 
on the brane and a five dimensional Gauss-Bonnet term in the bulk are present. In 
this realization, the Gauss-Bonnet term and the mass parameter in the bulk play a 
crucial role. 



I. INTRODUCTION 

The existence of dark energy is one of the most significant cosmological discoveries over 

n 

the last decade However, the nature of this dark energy remains a mystery. Various 
models of dark energy have been proposed, such as a small positive cosmological constant, 
quintessence, k-essence, phantom, holographic dark energy, etc., see Q] for recent reviews 
with fairly complete lists of references of different dark energy models. A cosmological 
constant is a simple candidate for dark energy. However, following the more accurate data 
a more dramatic result appears: the recent analysis of the type la supernovas data indicates 
that the time varying dark energy gives a better fit than a cosmological constant, and in 
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particular, the equation of state parameter w (defined as the ratio of pressure to energy 
density) crosses — 1 at z — 0.2 from above to below [3], where w = — 1 is the equation of 
state for the cosmological constant. The dark energy with w < — 1 is called phantom dark 
energy jj], for which all energy conditions are violated. To obtain w < —1, scalar field 
with a negative kinetic term, may be a simplest realization [a]. However, the equation of 
state of phantom scalar field is always less than —1 and can not cross —1. Also it has been 
shown that the equation of state cannot cross —1 in the k-essence model of dark energy under 
some reasonable assumptions f|. Some dark energy models which contain a negative- kinetic 
scalar field and a normal scalar field have been considered in [3]; in these models crossing 
the border w — — 1 can be realized. More recently it has been found that crossing —1 within 
a single scalar field model is possible, if its action contains higher derivative terms [q] (see 
also [i) . 

In this paper we suggest a possibility with effective equation of state crossing —1 in brane 
world scenario. The brane world scenario is now one of the most important ideas in high 
energy physics. In this scenario, the standard model particles are confined on the 3-brane, 
while the gravitation can propagate in the whole spacetime. As for cosmology in the brane 
world scenario, many works have been done over the last several years; for a review, see 
[lol ] and references therein. Brane world models admit a wider range of possibilities for 



dark energy 11 1. For example, the so-called DGP model has a late-time self-acceleration 
solution [l2|. On the other hand, higher derivative curvature terms, for example R 2 in 
Euler density, naturally appear in many occasions, such as in the quantum field theory in 
curved space Generally speaking, resulting equations of motion of such terms contain 
more than second derivatives of metric and the theory is plagued by ghosts. However 
there exists a combination of quadratic terms, called Gauss-Bonnet term, whose equation 
of motion has no more than second derivatives of metric and the theory is free of ghosts 
Q|. Another important property of Gauss-Bonnet term is that the Lagrangian is a pure 
divergence in four or less dimensions. So only in more than four dimensional theories, the 
Gauss-Bonnet combination has dynamically physical meanings. In addition , th e Gauss- 
Bonnet term naturally appears in low energy effective action of string theory |l5l . Ilq . The 
study on the effects of the Gauss-Bonnet term in the bulk of brane world model is therefore 
well motivated. 

The gravitational action can be generalized in various ways in brane world scenario. A 
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four dimensional scalar curvature term on the brane is an important one except a Gauss- 
Bonnet term in the bulk. This induced gravity correction arises because the localized matter 
fields on the brane, which couple to bulk gravitons, can generate via quantum loops a local- 
ized four-dimensional world- volume kinetic term for gravitons Q|. The combined effect of 
these curvature corrections to the action can remove the infinite-density big bang singularity, 
although the Gauss-Bonnet correction, which is expected to dominate at early times, on its 
own does not remove the infinite-density singularity, while the induced gravity correction 



on its own mostly affects the late-time evolution [171 ] . We find that the combining effect of 
the Gauss-Bonnet term in the bulk and the induced gravity term on the brane behaves as 
dark energy on the brane, and the effective equation of state can cross the phantom divide 
w = —1. 

Here we would like to mention that a dark energy model where standard gravity with 
dilaton scalar field containing an additional scalar-dependent coupling with Gauss-Bonnet 
invent has been investigated recently in Q; The effect the higher derivative terms on 
the big rip of phantom cosmology has been discussed in ll.9j. 

II. THE MODEL 

Let us start with the action 

S — "S'bulk + "Sbrane; ("0 

where 

Sbuik = Yq^g / rf5a V det((5) ^ ( (5)jR - 2As + aLcB ) ' (2) 



L GB =( 5 ) R* - 4 ®R AB ^R AB +( 5 > R ABGD {b) R ABCD , (3) 

and 

Sbrane = j d'xsjdet^g) (® R - 2A 4 ) . (4) 

Here S'bulk is the action of the bulk, Sbrane is the action of the brane, S is the total action, 
stands for the five dimensional Newton constant in the bulk, represents the four 
dimensional Newton constant on the brane, A 5 denotes cosmological constant in the bulk, 
A4 is cosmological constant on the brane, Lgb represents the Gauss-Bonnet term, and a is 
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a constant with dimension [length] 2 , is the induced metric on the brane 



(5) 



where n is the normal vector of the brane. We choose the coordinate of extra dimension y 
such that the brane stands at y — 0. Denoting 



( 4 )G" 



(6) 



one can see that r c has dimension of [length]. 

Assuming there is a mirror symmetry in the bulk, we have the Friedmann equation on 



the brane 



17(, see also 



— \l + -aH 2 + — + - 



H 2 + --U 



H 2 + 



k 8tt^G 



(P + A) 



(7) 



a* 3 

where H = a/a is the Hubble parameter, a is the scale factor, k is the spatial curvature of 
the three dimensional maximally symmetric space in the FRW metric on the brane, and 

A 4 



A 



8vr( 4 )G" 



(8) 



4a ^ 1 



;a 5 2m ( 5 )g 

4a ( + 



(9) 

6 a 4 ' 

Here M is a constant, standing for the mass of bulk black hole. Note that when one takes the 
po slt ,ve sign, the above e qU a tl on ca„ oe r ed u ced to the case o tg e n e rallZ e d DGP m ode, Q(M 
the pure DGP model, A4 = A5 = 0). In the latter case, we have 



4(8tt( 4 )G) 2 ( P + A)-2A 5 



tt2 k 8<k^G. x , 2 , 1 

H + a 2 - — ^ + ^^ 

Its various limits and inflation models have been studied in 



12 24M( 5 )G 



-.1/2 



. (10) 



2l| . Note that the branch with 
minus sign in (jHJ) cannot be reduced to the case of the generalized DGP model as a — > 0. 
Therefore we only consider the branch with the plus sign in (jOJ in what follows. 
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Introduce the following new variables and parameters, 

X =JP + MP = 1P~ fc ° ( + Z) ' 
8tt 

u = OTj2 (p + A) = Q m0 (l + zf + Q x , 
6H 

m = g«#0> 



n 



TJ2 2' 



l rrrr _ 2 1 / , f, 4^ 8aM( 5 )G\ 
=3^(- 1 + V 1 + — + ^?— J 

= 3^ (-1 + V / 1+^A 5 +^M0(1 + ^) 4 ) , (11) 

where ifo is the present value of Hubble parameter, z denotes the red shift with definition 
z = a Q /a — 1, and we have assumed that there is only pressureless dust in the universe. In 
addition, we have used the following notations 

o _ k o fa(4)(? fao o 87t(4)G X o - 3As o 3M(5)G 
ih °-~^Hl> ~ 3 Hfl' H$> ilA *-8Hl> ihl °-^Hf (12) 

With these new variables and parameters, (jZJ can be rewritten as 

4n(x - 2y)[l + m(x + y)} 2 = (x-u) 2 . (13) 

This is a cubic equation of the variable x. According to algebraic theory it has 3 roots. One 
can explicitly write down three roots. But they are too lengthy and complicated to present 
here. Instead we only express those three roots formally in the order given in Mathematica 

x\ = x±(y, u\m, n), 
%2 = x 2 (y, u\m, n), 

a?3 = x 3 (y,u\m,n), (14) 

where y and u are two variables, m and n stand for two parameters. The root on x of 
the equation ()13j) gives us the modified Friedmann equation on the Gauss-Bonnet brane 
world with induced gravity. From the solutions given in ()14|) . this model seems to have 
three branches. In addition, note that all parameters introduced in (jll)) and (|12jl are not 
independent of each other. According to the Friedmann equation p4jl . when all variables are 
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taken current values, for example, z — 0, the Friedmann equation will give us a constraint 
on those parameters 

1 = f(^w, ^mo, Qmo, ^a 5 , ^a, m, n). (15) 

Now we turn to the evolution of the universe in the brane world scenario, and pay 
particular attention to the acceleration characteristic of expansion at low red shift. To some 
issues on the dark energy, one of the most significant parameters from the viewpoint of 
observations is the equation of state, w = pressure/ energy density, of the dark energy. In 
general relativity the dark energy has to be introduced in order to coincide with observation 
data. In our model, the effect of the bulk Gauss-Bonnet term and induced gravity term on 
the brane may play the role of dark energy in general relativity. Therefore, in our model, 
the accelerated expansion of the universe is due to the combined effect of the brane world 
scenario and the Gauss-Bonnet term in the bulk. 

To explain observed accelerated expansion, we calculate the equation of state w of the 
effective "dark energy" caused by the bulk Gauss-Bonnet term and induced gravity term by 
comparing the modified Friedmann equation in the brane world scenario and the standard 
Friedmann equation in general relativity, because all observed features of dark energy are 
"derived" in general relativity. Note that the Friedmann equation in the four dimensional 
general relativity can be written as 

H 2 + - = ——(p + Pde ), (16) 

(2 o 

where the first term of RHS of the above equation represents the dust matter and the second 
term stands for the dark energy. Suppose that dark energy itself satisfies the continuity 
equation 

<h>de + 3H(p de +p de ) = 0, (17) 



dt 

we can then express the equation of state for the dark energy as 

p de ad\np de 

w de = — = -1 - -— . (18) 

Pde o da 

Note that we can rewrite the Friedmann equation (JUjl in the form of (JTfij) as 

tt2 8n^G f rr2 , x 8u^G \ 8ir^G . ^ , . 

xHl = — — p+ ( H*x(y,u\m,n) — p\ (p + F), (19) 
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where p is the energy density of dust matter on the brane and the term 

F = ( ^f^ G x{y,u\m, n) - p 

corresponds to pd e in ([TE]). In addition, here x in (TH?]) denotes one of three solutions given 
in (|14|). Now we regard F in (|19|) as the effective dark energy density, then we have the 
equation of state for the effective dark energy 

, adlnF _ 1 dlnF 

W ef f = -1 : = -1 + -— — -. 20 

/; 3 da 3d\n(l + z) K ' 

The type la supernovas observation data indicate that the equation of state of dark energy 
undergoes a transition from w > — 1 to w < — 1 at z = 0.2 [3J]. As pointed out above, it is 
quite difficult to realize such a transition in usual field theory models. Here we will see that 
it can be realized in our Gauss-Bonnet brane world scenario with induced gravity, in which 
the Gauss-Bonnet term and the mass parameter M in bulk play an important role. 

Before proceeding, let us first discuss the case without the Gauss-Bonnet term, namely 
a = 0. In this case, as can be seen from fllOL the effective dark energy density is 



2 _ j 



' / us„\2, . . 12 24M^G ll/2 ' 
A(8n^G) ( p + A)-2A 5 + — 



r£ a 4 



:2D 



SttWG [r c 2 v/3r 

When the bulk mass parameter M vanishes, namely, M = 0, this form of effective dark 
energy has been discussed in For the branch with plus sign in (|21|) ( called branch 2 
in the equation of state for the effective dark energy is always larger than —1, and 

approaches to — 1 as z — > — 1. On the other hand, for the branch with minus sign in (|21j) ( 



called branch 1 in lip, there is a divergence on the curve of w e ff versus the red shift z at 
some point z , and in the region z G (—1, z ), w e ff monotonically goes from — oo at z = z 
to — 1 at z — —1, while in the region z G (zq, oo), it monotonically goes from +oo at z to 
some value larger than — 1 at z = oo. Therefore in the branch 1, the equation of state for 
the effective dark energy can be less than —1, but it can never cross the phantom divide — 1. 
In addition, let us remind that the divergence of the equation of state at some point z = zq, 
whose value depends on the model parameters, is harmless. This divergence is caused by 
F = there, except for the equation of state w e ff, other physical quantities like Hubble 
parameter, deceleration factor, energy density of the effective dark energy, etc., behave well. 

When the bulk mass parameter M does not vanish, namely M^O, one may ask whether 
this situation will be changed. The answer is no. In this case, the equations of state for two 
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branches are 

+ . 3Ia-8J 



w + = -1 + 



QB(a 2 + B) 



31a — 8 J , . 

ur = - 1 - m^wy (22) 

where w + and w~ stands for the " + " and " — " branch in (J2TJ), respectively, and 



I = Pma(4, B = Va 4 + al - 2J>, J = Ar~ 2 (5) G. 

Note that w + diverges at B = 0, and w~ diverges at a 2 = B or B = 0. We further find from 
that if B ^ 0, we have w + = —1 and w~ = —1 when 

8J 



However, we can show that 

dw + 



da 



a = — . 
3/ 



dw 

> 0, 



da 



> 0. 



Both cases are not consistent with the observation fact: the equation of state of dark energy 
transits at z — 0.2 from w > —1 for z > 0.2 to w < —1 for z < 0.2 0. 

Next we turn to the case with q ^ 0, namely the case with a Gauss-Bonnet term in the 
bulk. For simplicity, we set k = and A4 = A5 = in what follows. To see the properties 
of three branches in (|14jl. let us first look at the asymptotical behavior of them as z — *■ — 1. 
In this case, the equation (fT3*j) reduces to 



Anx{l + mxf = x 2 , (23) 



which has three roots 



Xi = 0, 

1 — 8nm + y/1 — 16nm 



x 2 



inm 2 



1 — 8nm — \/l — lQnm , A . 

^3 = ^ 5 • 24 

We find that when nm > 1/16, three solutions in (fT4")) go to the limits in ()24j) in the following 

manners 

xx -> xi, x 2 -> x 2 , x 3 -> x 3 ; (25) 
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when nm < 1/16, 

Xi x 2 , x 2 -> xx, x 3 -> x 3 ; (26) 

and when nm = 1/16, 

x\ — > 1/m, X2 — > 0, S3 — > 1/m. (27) 

From the above, we can see that (1) when nm > 1/16, only does the solution x\ make 
sense, while x 2 and x% are not physical solutions since x 2 and x 3 are imaginary; (2) when 
nm < 1/16, three roots Xi, x 2 and 23 in (j24j) look reasonable since they are all real. But we 
note that when m — > 0, the branches x 2 and X3 reduce to two branches in the DGP model, 
while x\ — » 00. So is a new branch in our model and has no corresponding one in the 
DGP model. The two branches in the DGP model also can be obtained directly from (JHJ): 



lim x 2 = u + 2n — 2\ / n 2 + 



rax, 



m— >0 



lim £3 = u + 2n + 2\J n 2 + nu. (28) 
We now consider the case M = in our model. Assuming Q m0 = 0.28, which is indicated 



by WMAP data 



22|, we find that for any values of m and n, the condition 

x\ z=0 = 1, (29) 



always cannot be satisfied for the branches xi and x 2 in (JHJ), while X3 always becomes 
imaginary at finite red shift. Therefore the case with a vanishing bulk mass, M — 0, cannot 
satisfy our requirement. In Fig. 1 we plot the x for the branch x± with respect to the 
parameters m and n at red shift z — 0. 
The case with M^0: 

(i) For the branch Xi in (pPfj) , we find that it is possible to satisfy the condition Xi| z= o = 1 
with f2 m o = 0.28. From the viewpoint of physics consistency we may think from (fTTj) that 
m and n should be of order one or less. And from (fTK|) one can obtain 

Rmo = tt MQ (m,n), (30) 

with condition xi| 2= o = 1 and Q m o = 0.28. In Fig. 2 we plot Qmo as a function of m and n, 
which indicates that Qmo is positive and is of order one, as expected. 

In Fig. 3 we show the equation of state for the effective dark energy when we take 
m = 1.036 and n = 0.04917. In this case, from the constraint equation ()15|) . one has 
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FIG. 1: xi| z= o with M = and ft, m o = 0.28 versus the parameters m and n. It is obvious that 
x\ (m,n) is not a continuous function. The discontinuity happens at nm = 1/16. In the left part 
of the surface x(m,n) S> 1, the right part x(m,n) <C 1. There does not exist a point on which 



Therefore the equation of state for the effective dark energy can indeed cross the phantom 
divide w — — 1 near z ~ 0. 

Some remarks are in order. First when nm > 1/16, one can also realize crossing w — — 1 
for the equation of state in this model. However, in this case, one may find that the 
effective dark energy density F in (|19jl is always negative, which cannot drive the universe 
to accelerated expand. Therefore the stability condition nm < 1/16 has to be satisfied. 
Second, as the case of generalized DGP model, the equation of state w e ff in this model has 
also a divergence at the some red shift point, which depends on the model parameters. The 
divergence is harmless and it is also caused by F = at that point. To see the harmlessness, 
In Fig. 4 and 5 we respectively plot the Hubble parameter and the fraction energy density 
of dust matter in the case of m — 1.036 and n = 0.04917. From Fig. 4 we can see that 
when z < 0, the Hubble parameter is an increasing function of red shift, which shows the 
phantom behavior of dark energy. The universe finally approaches to a de Sitter space. 



xi = 1. 



VLmo = 2.08. From the figure we see that w e // < — 1 at 2; = and 



d,W e ff 



< 0. 



dz 



2=0 
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FIG. 2: £Imo versus m and n, with J7 m o = 0.28 and k = 0. 
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FIG. 3: The equation of state w e ff with respect to the red shift 1 + z, with {l m o = 0.28 and 
n A = 2.08. 

To further show the harmlessness of the divergence in the equation of state w e ff, we 
would like to plot in Fig. 6 the deceleration parameter q, which is denned by 

q=~ = ±[Q m + (l- n m ){l + 3w eff )}. (31) 

The transition of the universe from deceleration to acceleration phase happens at z — 0.8 
in this parameterized model. 

In the above demonstration, we take Q m0 = 0.28. In Fig. 7,8 and 9, we further plot the 
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FIG. 4: # 2 /#o versus 1 + z ' with °mO = 0.28 and ft A = 2.08. 
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FIG. 5: Q m versus 1 + z, with f2 m o = 0.28 and fi^ = 2.08. 

corresponding quantities for the case with Q m0 = 0.05, m = 1.17 and n = 0.0511. 

(ii) For the branches X2 and 23 in (I14j) . we find that they always become imaginary at some 
finite red shift. It can be understood as follows. These two branches go to corresponding 
ones of the generalized DGP model as m — > 0, respectively. For the generalized DGP model, 
we can see from (J2T|) that there is a square root in the effective dark energy density and the 
square root will always become imaginary as a decreases if M 7^ 0. 
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FIG. 6: q versus l + z, with O m0 = 0.28 and U A = 2.08. 
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FIG. 7: w e ff versus l + z, with O m0 = 0.05 and = 0.407. 

III. CONCLUSION 

To summarize, in this paper we have shown that in the Gauss-Bonnet brane world with 
induced gravity, where a Gauss-Bonnet term in the bulk and a scalar curvature term on the 
brane appear, it is possible to realize crossing the phantom divide w = — 1 for the equation 
of state of the effective "dark energy", without introducing any phantom component. The 
dark energy model can be consistent with astronomical observation data. In realizing the 
transition from w > — 1 to w < - 1, the Gauss-Bonnet term and the bulk mass parameter 
play a crucial role. Finally we point out that the qualitative conclusion will not be changed 
if we include non-vanishing cosmological constants on the brane and in the bulk, namely, 
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FIG. 8: H 2 /Hl versus 1 + z, with VL m0 = 0.05 and £l A = 0.407. 

Omega_m 
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FIG. 9: n m versus 1 + z, with O m0 = 0.05 and = 0.407. 

A 4 7^ and A 5 ^ 0. In addition, it would be of great interest to fix the parameters in this 
dark energy model, according to current type la supernovas data. This issue is currently 
under investigation. 
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